APS/123-QED 



Energy-dependent Lorentz covariant parameterization 
of the NN interaction between 50 and 200 MeV 

Z. P. Li 1 , G. C. Hillhouse 2 ' 1 , and J. Meng 1 - 2 - 3 ' 4 
1 School of Physics, Peking University, Beijing 100871 
2 Department of Physics, University of Stellenbosch, Stellenbosch, South Africa 
^Institute of Theoretical Physics, Chinese Academy of Science, Beijing 100080 and 

4 Center of Theoretical Nuclear Physics, 
National Laboratory of Heavy Ion Accelerator, Lanzhou 730000 
(Dated: February 2, 2008) 

Abstract 

For laboratory kinetic energies between 50 and 200 MeV, we focus on generating an energy- 
dependent Lorentz covariant parameterization of the on-shell nucleon-nucleon (NN) scattering am- 
plitudes in terms of a number of Yukawa-type meson exchanges in first-order Born approximation. 
This parameterization provides a good description of NN scattering observables in the energy range 
of interest, and can also be successfully extrapolated to energies between 40 and 300 MeV. 
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I. INTRODUCTION 



We present an energy-dependent Lorentz covariant parameterization of the on-shell 
nucleon-nucleon (NN) scattering matrix at incident laboratory kinetic energies ranging from 
50 to 200 MeV. In particular, we focus on the relativistic Horowitz-Love- Franey (HLF) model 
which parameterizes the NN scattering amplitudes as a number of Yukawa-type 
meson-exchange terms. 

The original HLF model parameterized pp and pn scattering amplitudes at discrete en- 
ergies of 135, 200 and 400 MeV A drawback of the initial representation was that 
separate fits were performed at each energy resulting in non-systematic and unphysical 
trends in the various model parameters as a function of energy. This was also undesirable in 
the sense that it hindered meaningful comparisons of the off-shell properties of the NN scat- 
tering matrix at different energies. This problem was subsequently addressed by Maxwell 
who developed energy- dependent parameterizations of the NN scattering amplitudes for two 
energy ranges, namely 200 - 500 MeV jsj and 500 - 800 MeV Q|. Since several experiments 
at future radioactive ion beam facilities will be conducted at nucleon energies lower than 
200 MeV, it is essential to extend the HLF model to lower energies so as to generate reliable 
input for nuclear reaction models. 

An attractive feature of HLF model is the existence of a simple relationship - lacking in 
conventional nonrelativistic models - between the Lorentz invariant NN scattering ampli- 
tudes and mesons mediating the interaction. Moreover, this model has a physical basis in 
the one-boson-exchange (OBE) picture, since the values of the real meson-nucleon coupling 
constants (at energies higher than 200 MeV) are similar to those arising from more sophis- 
ticated OBE models [1]. In addition to being employed as the basic NN interaction driving 
the dominant reaction mechanism in several relativistic scattering models of nucleon- induced 
reactions, the HLF model has also been applied to extract relativistic microscopic optical 
potentials (3, O, 6, Q] for studying elastic proton scattering from stable nuclei, as well as 
for generating scattering wave functions for evaluating transition matrix elements in various 
relativistic distorted wave models {5, f|. More specifically, the latter microscopic optical po- 
tentials are generated by folding the HLF t-matrix with the relevant relativistic mean field 
Lorentz densities via the so-called tp-approximation. Indeed, one of the triumphs of this 
folding procedure is that the microscopic potentials, for energies from 200 to 400 MeV, are 
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virtually identical to the corresponding global phenomenological Dirac optical potentials, 
which have been shown to provide excellent quantitative predictions of scattering observ- 
ables for elastic proton scattering from spin-zero stable nuclei ranging from 12 C to 208 Pb and 
for incident energies between 20 and 1040 MeV 8|. 

One of the frontier areas of research in nuclear physics is to understand the properties 
of nuclei far from the beta-stability line, the so-called exotic or unstable nuclei. These 
studies provide insight into the nuclear processes that underlie the evolution of the stars 
and the origin of the elements in the cosmos. Indeed, a number of radioactive ion beam 
facilities are currently under construction to study short-lived rare isotopes. In particular, 
there are plans to study exclusive proton-induced proton-knockout reactions (using inverse 
kinematics) at facilities such as RIKEN and GSI for both neutron- and proton-rich nuclei at 
energies lower than 200 MeV per nucleon. For the eventual interpretation of these data it 
is essential to have reliable optical potentials. Existing global Dirac optical potentials have 
been constrained to reproduce elastic proton scattering from stable nuclei jsj], and hence 
there is no reason to believe that these potentials can be reliably extrapolated for studies 
of unstable nuclei. On the other hand, one can readily extend the above-mentioned folding 
procedure to calculate scattering potentials for proton-induced reactions on exotic nuclei. 
Moreover, the successful application of the relativistic tp approximation to describe elastic 
proton scattering from stable nuclei, gives one confidence in extending this approach to study 
proton scattering on exotic nuclei. Two basic ingredients underly the realization of these 
folding potentials, namely a suitable analytical representation for the NN interaction in the 
energy range of interest, as well as an appropriate relativistic model of nuclear structure 
Isjl. Several relativistic mean field nuclear structure models are currently being developed 

n n 

for studying the structure of unstable nuclei [9|, [10J. In this paper, however, we focus 
on developing a HLF-type parameterization of the NN scattering amplitudes at energies 
lower than 200 MeV. The generation of the corresponding microscopic optical potentials will 
be considered in a future paper, in which Pauli blocking corrections and nuclear medium 
modifications to the NN interaction are expected to be significant and will be taken into 
account 12|. 



The HLF model is described in Sec. (jlT]). The fitting procedure and corresponding results 
are presented in Sees, f illip and (II VI) . respectively. 
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II. THE RELATIVISTIC HOROWITZ-LOVE-FRANEY MODEL 



The most general parametrization of the nonrelativistic on-shell elastic NN scattering ma- 
trix - consistent with rotational, parity, time-reversal and isospin invariance - is completely 
determined by five complex scattering amplitudes. Depending on the specific application of 
interest, the NN scattering matrix can be recast into several distinct, but equivalent, five- 
term representations, one of which is the so-called Bystricky parameterization [3] , namely: 

(21k)- 1 f cm (k, k') = ^{(a + b)h ® I 2 + (a - b)(a x <g> I 2 ) • h(h ® a 2 ) ■ ft + (c + d) 
x (cti ® h) ■ m{h ®a 2 )-m+(c- d){ffi ® h) ■ £{h ® o 2 ) • I 
+ e(ai ® I 2 + A ® a 2 ) ■ n} (1) 

where ® denotes a kronecker product, <?i and a 2 refer to the usual Pauli spin matrices 
associated with the projectile and target nucleons, 1 and 2 respectively, I\ and I 2 are the 
respective 2x2 unit matrices, and the basis vectors £, rh and n are defined as: 

rn = fi = , (2) 



\k' + k\ ' \k' - k\ ' \k x fc'| ' 

where and k! are the initial and final momenta of the interacting nucleons in the NN centre- 
of-mass frame. The scattering matrix is normalized such that the polarized differential cross 
section for free NN scattering is given by 

(^)pol = \(xl 1 Xi 2 \fcm\Xs 1 Xs 2 )\ 2 (3) 

where the x's represent the usual Pauli spinors. 

For relativistic applications, the preceding nonrelativistic phenomenology, expressed by 
Eq. ([1]), can be cast into many different relativistic Lorentz invariant forms of which the 
on-shell NN matrix elements are identical to those of the nonrelativistic scattering matrix. 
A convenient choice of the relativistic NN scattering matrix, commonly referred to as the 

n 

IAl representation and originally proposed by McNeil, Ray and Wallace [L4J, is given by: 

5 

i NN = Y, F L(^u)\ i l ) -\ ( ? (4) 

L=l 

where s, t and u quantities denote the usual Mandelstam variables, s = (ki + k 2 ) 2 , t = 
(ki — k[) 2 , where u = (k\ — k' 2 ) 2 , and k\ = (E,k) and k 2 = (E,—k) define the incident 
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4-momenta of the projectile and target nucleons respectively in the center-of-mass frame, 
with similar definitions holding for the outgoing (primed) nucleons. The L ,s index the five 
Dirac matrices listed in Table HJ and the dot product implies contraction of the Lorentz 
indices. 



TABLE I: Dirac matrix types parameterizing the free NN amplitudes. 



L 


Lorentz type 


Dirac matrices 


1 


Scalar (S) 


^4x4 


2 


Vector (V) 


%< 


3 


Pseudoscalar (P) 


75 


4 


Axial-vector (A) 


757^ 


5 


Tensor (T) 





The transformation between the nonrelativistic on-shell NN a, b, c, d, e amplitudes and 
the relativistic F s , Fy, F P , Fa, F t amplitudes is readily established by equating the matrix 
elements of the nonrelativistic scattering matrix, given by Eq. (HI), to the matrix elements 

nnnri 

of the relativistic scattering matrix, given by Eq. gl) [1|, 12, la |4fl, namely: 

(2iky 1 xl l xl f cm X Sl Xs 2 = u(k[,s' 1 )u(k' 2 ,s' 2 )i NN u(k 1 ,s 1 )u(k2,s 2 ) , 



(5) 



where the Dirac spinors u are normalized such that uu = 1 . The corresponding 5x5 matrix 
(9 5x5 linking the two representations is given by 

fa\ f Fs \ 



b 

c 
d 



ikO, 



5x5 



Fy 
F P 
F A 
\F T J 



(6) 



where 



O 



5x5 



/ 


cos# 


cos 6 


cos 9 


— cos 9 


^sin^ ^ 


( 


a 


ea + p/3 





pa + e/3 


-2(3 


\ 




1 


-1 


1 


1 







— e7 


-7 


-VI 


-5 


2(p + e)6 






-1 


1 


1 


1 







a 


(p + e)a 





(e — p)a 


2a 






1 


1 


-1 


1 







e7 


7 


-vi 


-7 


2l 




V 


—i sin 9 


— % sin 9 


—i sin 9 


i sin 9 


— 4i cos# J 


V 


— e 


— e 





e 


-2e 


/ 



and 



01 = 


2^ 

cos — 
2 


p = 


1 + snr - 

2 


7 = 


sm - 
2 


5 = 


2# 

1 + cos - 
2 


e = 


e e e 

— sm - cos - 
M 2 2 


E = 


(k 2 + M 2 ) 1 / 2 


e = 


£ 2 
M2 


P = 


k 2 



with being the NN centre-of-mass scattering angle and M the free nucleon mass. The 
HLF model is described in detail in Refs. Here we briefly allude to the main 

aspects relevant to this paper. Essentially this model parameterizes the complex relativistic 
amplitudes Fi(s,t,u) in terms of a set of N = 10 meson exchanges [see Table Hi] in first- 
order Born approximation, such that both direct and exchange NN (tree-level) diagrams are 
considered separately, that is: 

iM 2 



F L (s,t,u) = ^[FE(s,t) + FX(s,u)} , (7) 



where 



N 



F?(s,t) = ■ ^} Ti r(£, l9 I) (8) 



i=l 



N 



f*( s ,u) = (-i) T ™Y, c mMn-^) Ti nE,\Q\) ■ (9) 



1=1 



Here = (0,1) denotes the isospin of the i meson, TjviV refers to the total isospin of the 



:ix 



two-nucleon system, Cl(i),l is the Fierz matrix [1|, [15J, and 

f(E,x) = r R {x)-if i I (x) (10) 

with 
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fKx) = ^— 2 a + f 2 r\ (12) 
x z + ml Af 

where x represents either the direct three-momentum transfer \q\ or the exchange- 
momentum transfer \Q\, namely 

q = k — k', 

Q = k + k' . (13) 
The isospin matrix elements in Eqs. (jHJ) and (jHJ) yield 

„ 1 for the exchange of T, = (isoscalar) mesons , 

(ri ■ r 2 ) T ^ = J } (14) 

1 for the exchange of T{ = 1 (isovector) mesons 

for the Tnn = 1 amplitudes and 

T 1 for the exchange of T« = (isoscalar) mesons 

(n ■ t 2 ) 1 = < } (15) 

—3 for the exchange of Tj = 1 (isovector) mesons 

for the Tjvjv = amplitudes. 

Note that the mesons in Table HT1 represent different Lorentz types (S, V, P, A, T) with 
an isospin of either or 1. The coupling constants are complex with real and imaginary 
parts, gl and gf, respectively. The imaginary couplings represent a purely phenomenological 
means of parameterizing the imaginary amplitudes. The meson propagators are 

1 

x 2 + m 2 ' 

and the following monopole form factors are assumed for the meson-nucleon vertices, 

1 



(16) 



1 + 



(17) 



with separate masses (m^, rrii) and cutoff parameters (Aj, Aj) for the real and imaginary 
parts of the amplitudes, respectively. 



III. FITTING PROCEDURE 



Using some initial guess (see later) for the HLF model parameters, the relativistic F$, 
F v , F P , Fa, F t amplitudes are determined [from Eq. ((JJ)], and these in turn are converted 
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to the nonrelativistic a, b, c, d and e amplitudes- via Eq. - and subsequently compared 
to the SP05 (Spring 2005) empirical amplitudes to establish the goodness of the fit. More 
specifically, the HLF parameters are varied to simultaneously fit the amplitudes a, b, c, d, e 
for both Tjvtv = (0, 1) at laboratory kinetic energies between 50 and 200 MeV (in 25 MeV 
steps) and centre-of-mass scattering angles between 5 and 175 degrees (in 5 degree steps). 
Empirical values of the a, b, c, d and e amplitudes were extracted via the on-line Scattering 
Analysis Interactive Dial-in (SAID) facility. On a more technical note, the NN amplitudes 
were obtained via ssh call to the SAID facility, gwdac.phys.gwu.edu, with user id said (no 
password), in which there are 4 choices for different types of NN scattering, namely pp, np, 
npl, npO. For this paper we selected the npl and npO options for generating the T^n — 1 
and T/vat = amplitudes, EXCLUDING Coulomb corrections, respectively. After choosing 
the isospin of interest, one is confronted with 6 different choices for different nonrelativistic 
amplitude types: 1 (VPI-H), 2 ( Wolfenstein) , 3 (Bystricky), 4 (Helicity), 5 (Transversity), 
or 6 (Transverse planar). We chose option 2 (Wolfenstein) which actually represents the so- 
called Hoshizaki amplitudes [3]. The latter dH, rnH, gH,hn and ch Hoshizaki amplitudes 
are subsequently converted to the required Bystricky a, b, c, d and e amplitudes [in Eq. (CQ)] 
using the transformation [l^ : 

a = aH + rriH, b — aH — rnu, c = 2gn, d=—2hn, e = 2ch ■ (18) 

Separate fits were performed to the real and imaginary parts of the amplitudes. In order for 
the 5 empirical amplitudes to be weighted equally in the fit, we minimized the value of x 2 
defined by 

2 V v (-^empirical 

data X empirical/ 

employing a Levenberg-Marquadt method. Here (^empirical) represents an angle-averaged 
value. For a specific amplitude, say a, and T lab the angle-averaged value is determined by: 

175° 

5-^ ^empirically) 

( Q empirical) = (20) 
-< * ang 

where iV an g = 35 represents the number of angles fitted per energy. The total number of 
the data points per real or imaginary fit is 2450: 7 energies x 35 angles x 5 amplitudes x 
2 isospins. 
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FIG. 1: (Color online) The experimental SP05 pseudoscalar amplitudes Tnn = 1 F p , represented 
by the solid squares, as a function of the laboratory kinetic energy Ti a b for a small centre-of-mass 
scattering angle of 5 degrees. The solid curve is an exponential function of the form given by 
Eq. (EH). 



We now discuss the general considerations for extracting the optimal HLF coupling 
strengths, meson masses and cutoff parameters. The energy dependence of the coupling 
constants was established as follows. The dominant contribution to the pseudoscalar Fp 
amplitude is determined almost completely by pion exchange to the direct term [l], and 
thus according to Eqs. Q - (fT2l) the amplitude F P at small centre-of-mass scattering an- 
gles, and hence small values of \q\, is approximately proportional to g\. Consequently the 
relationship between F v and Tl a b (the laboratory kinetic energy) is essentially the same as 
the functional dependence of g\ on Ti ab . In Fig. [fl we plot the F P amplitude (for T NN = 1) 
as a function of Tj a b for a small centre-of-mass scattering angle of 5 degrees. A clear expo- 
nential energy-dependence is observed. Guided by the latter trend, we chose the following 
exponential energy- dependence for all (both real and imaginary) coupling constants, namely 



g\E) = g 2 [l + a g {e 



a T T r( 



1)] 



(21) 



where 

T rcl ee (22) 

with To = 200 MeV, T rc i is positive in the 50 to 200 MeV energy range of interest, and 
<7o; a g and fl r are dimensionless parameters extracted by fitting to the relevant data. For 
the real meson masses we chose the experimentally measured values 16[ for the ir, u and 
p mesons, namely 138, 782, and 770 MeV, respectively. For the other real and imaginary 
meson masses, as well as the real and imaginary cutoff energies A, and Aj, we chose the 
starting the values to be the same as those published in Table I of Ref. 

Firstly, the meson masses, coupling constants and cutoff parameters were varied slightly 
to obtain the best fit to the 200 MeV data alone. Similar to Ref. jsj], the values of Aj and 
Aj were restricted to always exceed the respective meson masses. Thereafter, the meson 
masses and cutoff parameters were kept fixed, and a g and or, which determine the coupling 
constants via Eq. ff2TT) . were varied to obtain the best fit to the total data set between 50 
and 200 MeV. Then all parameters (with the exception of the meson masses) were varied to 
fit the total data set. 



IV. RESULTS 

The fits to the real and imaginary amplitudes yield minimum \ 2 values of 10.203 and 
10.798, respectively: the relative \ 2 values can be obtained by dividing the latter values by 
2450. The fit to the real part of amplitude a for Inn = was found to be inferior compared 
to the real parts of the other amplitudes. Therefore, for the real part, the weight of amplitude 
a was multiplied by 200, and the new x 2 value is equal to 16.657. The imaginary parts of 
amplitudes b, c for Inn = 1 and a, b, c for Inn = showed systematic deviations in the low 
energy range. For the imaginary amplitudes, the parameter cjt of the isovector tensor meson 
was changed slightly, such that the systematic deviation nearly disappeared, and the fits to 
the other amplitudes remained satisfactory. Fitted values for the various HLF parameters 
are presented in Table [III 

To demonstrate the quality of the fits we compare, in Figs. (j2J) to (j4)), the fitted amplitudes 
a, b, c, d, and e, for both Tjvjv = (0, 1), to the corresponding SP05 empirical values at selected 
energies of 50, 125 and 200 MeV. Our HLF parameter set provides a satisfactory description 
of the empirical amplitudes at the latter energies. Although not displayed, this is also the 
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TABLE II: Real and imaginary HLF parameters. The masses m and cutoff parameters A are in 
MeV, whereas the other parameters are dimensionless. 



Real parameters 



Meson 


Isospin 


Coupling type 


m 


9l 


a 9 


ClT 


A 


a 





Scalar (S) 


600 


-8.379 


-5.581 x 10" 1 


2.364 x 10" 1 


965.23 


CO 





Vector (V) 


782 


1.014 x 10 1 


-2.219 x 10~ 2 


3.793 


1158.74 


to 





Tensor (T) 


550 


2.783 x 10- 1 


1.331 


3.074 


1955.59 


a 





Axial vector (A) 


500 


4.842 x lO" 1 


1.440 


2.847 


1577.53 


V 





Pseudoscalar (P) 


950 


1.089 x 10 1 


1.553 


1.956 


980.99 


s 


1 


Scalar (S) 


500 


6.233 x 10~ 3 


5.675 


3.521 


3000.00 


p 


1 


Vector (V) 


770 


-1.535 x lO" 1 


3.429 


1.825 x lO" 1 


3000.00 


h 


1 


Tensor (T) 


600 


-2.497 x 10" 1 


5.508 x lO" 1 


3.959 


1290.72 


a i 


1 


Axial vector (A) 


650 


-1.355 


2.480 x 10" 1 


3.623 


745.19 




1 


Pseudoscalar (P) 


138 


1.195 x 10 1 


-1.671 x 10" 1 


3.216 x 10" 1 


678.44 


Imaginary parameters 


Meson 


Isospin 


Coupling type 


m 


9l 


a g 


a-f 


A 


a 





Scalar (S) 


600 


-2.866 


7.722 x 10" 1 


1.503 


772.74 







Vector (V) 


700 


4.415 


7.094 x 10" 1 


1.528 


701.00 


to 





Tensor (T) 


750 


-9.114 x 10" 1 


1.217 


1.674 


1395.12 


a 





Axial vector (A) 


750 


-2.124 


1.187 


1.524 


1364.90 


V 





Pseudoscalar (P) 


1000 


1.411 x 10 1 


3.021 x 10~ 2 


4.970 


3000.00 


s 


1 


Scalar (S) 


650 


3.089 


9.873 x lO" 1 


8.292 x 10" 1 


771.14 


p 


1 


Vector (V) 


600 


-2.464 


7.620 x lO" 1 


1.176 


795.81 


h 


1 


Tensor (T) 


750 


6.447 x lO" 1 


3.316 x 10- 1 


2.977 


1741.21 


a i 


1 


Axial vector (A) 


1000 


2.328 


3.885 x lO" 1 


2.915 


1256.95 


7T 


1 


Pseudoscalar (P) 


500 


-5.788 


7.635 x 1Q- 1 


1.623 


1391.26 
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50 MeV 




e (deg) 

FIG. 2: (Color online) Comparison of the fitted real and imaginary (a,b,c,d,e) Tjvjv = 1 and 
Tnn = amplitudes (dashed curve), in units of fm, with the empirical SP05 amplitudes (solid 
curve) as a function of the NN centre-of-mass angle 9, in degrees, at a laboratory kinetic energy 
50 MeV. 

case at all other energies between 50 and 200 MeV. 

The quality of our fit is also judged by comparing pp and pn HLF-based and empirical 
scattering observables derived from the a, b, c, d, e amplitudes. In order to calculate the rele- 
vant pp and pn observables it is necessary to extract the corresponding pp and pn amplitudes 
via the following relation: 

F^pp) = F t (T = l) (23) 



Fi(pn) = ^mT = l) + F t (T = 0)], 
12 



(24) 



125 MeV 



— s 0.3 
0.0 

(Q 

-0.3 
0.8 

0.0 
0.0 

o -°- 4 

-0.8 
0.6 

T3 

0.3 
0.008 

<D 0.000 

-0.008 




T NN = 1 



I I I I I I I 




I I I I I I I 



I I I ■ I ■ I 




0.2 
0.0 
-0.2 

0.3 

0.0 

0.0 

-0.3 
0.0 

-0.1 
0.6 

0.0 



T NN = 1 



Imaginary 
I ■ I ■ I ■ I 



I ■ I ■ I ■ I 



I i I i I ■ I 




I i I i I i I 




1.5 



1.0 

2 

1 


2 
1 

1 



-1 

0.3 

0.0 
-0.3 




I ' I ■ I ■ I 



I i I i I ■ I 




I ■ I ■ I ■ I 




1.0 
0.5 
0.8 
0.4 

0.0 
0.8 

0.4 

0.0 
0.1 

0.0 

-0.1 
0.2 
0.0 
-0.2 




Imaginary T NN = 
I ■ I ■ I ■ I 




60 120 180 60 120 180 60 120 180 

e (deg) 



i ■ i ■ i 




60 120 180 



FIG. 3: (Color online) Comparison of the fitted real and imaginary (a,b,c,d,e) Tnn = 1 and 
Tnn = amplitudes (dashed curve), in units of fm, with the empirical SP05 amplitudes (solid 
curve) as a function of the NN centre-of-mass angle 9, in degrees, at a laboratory kinetic energy 
125 MeV. 

and then convert them to a,b,c,d,e amplitudes via Eq. (Q. For eventual comparison to 
experimental observables, one must add Coulomb corrections (which dominate at forward 
and backward scattering angles) to the a,b,c,d,e amplitudes for Tnn = 1 ; the Coulomb 
amplitudes were obtained by subtracting the SAID Tnn = 1 amplitudes from the SAID pp 
amplitudes. The NN scattering observables of interest are defined as follows [l]: 
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FIG. 4: (Color online) Comparison of the fitted real and imaginary (a, b,c, d, e) Tjvjv = 1 and 
Tnn = amplitudes (dashed curve), in units of fm, with the empirical SP05 amplitudes (solid 
curve) as a function of the NN centre-of-mass angle 9, in degrees, at a laboratory kinetic energy 
200 MeV. 
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A 



a 
P 
D 

yy 
R 

A 



■(|a| 2 + |6| 2 +|c| 2 + |rf| 2 + |e| 



Re(a*e)/cx 



1 

2 
A 

D, 



M 2 +|e| 2 )Ax 



(|a| 2 -|6p 



+ |rf| 2 + |e| 2 )/ C r 



$ $ 6 

Re(a*b) cos(a + -) + Re(c*d) cos(a — -) — lm(b*e) sin(a + —) 

6 6 6 

— Re(a*b) sin(a + -) + Re(c*d) sin(a — -) — lm(b*e) cos(a + —) 



1° 
1° 



where 



6 

a= 2 



'lab 



(25) 



with 6*i a b being the laboratory scattering angle, and a = da/dQ represents the unpolarized 
differential cross section. 

Before comparing empirical NN scattering observables to our optimal HLF predictions, 
we briefly comment on the parameter-sensitivity of the latter. We have established that the 
lower energy observables are more sensitive to variations in the values of g 2 [calculated from 
(?q, a g and in Eq. (I2ip ]. For example, for pp scattering at 50 MeV, variations of 1 % 
in the g 2 values translate to a maximum change of 6 % on the polarization (P), the effect 
being reduced for all other spin observables, and also modify the unpolarized differential 
cross section (da/dQ) by 15 %, where the latter corresponds to a 1 % change in the real 
isoscalar tensor coupling. On the other hand, all the pn observables exhibit less sensitivity. 

Results for the unpolarized differential cross section (da/dQ), the polarization (P), the 
depolarization (D), the tensor asymmetry (A TO )and triple scattering parameters A and R 
are shown in Fig. ([5]) to (J7J) for pp and pn at three different energies, namely 50, 125 and 
200 MeV. In general, our HLF parameterization provides an excellent description of the 
empirical scattering observables. In particular, the quality of our fits is just as good, if not 
better, than those presented in Refs. . Although this paper focussed on the 50 to 

200 MeV range, we have also confirmed that our parametrization provides a satisfactory 
description of scattering observables at a lower enery limit of 40 MeV and a higher energy 
limit of 300 MeV. The next phase of this project will be to test the predictive power of 
the relativistic tp folding procedure (for generating microscopic relativistic scalar and vector 
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optical potentials) for describing elastic proton scattering from nuclei at energies lower than 
200 MeV. Indeed, at these low energies, Pauli blocking corrections and nuclear medium 
modifications to the NN interaction will play a significant role (2), ll|, 12]. RIA predictions as 
well as the importance of the latter corrections for elastic proton scattering will be presented 
in a forthcoming paper. 
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FIG. 6: (Color online) pp and pn scattering observables versus the centre-of-mass angles 6 at an 
incident laboratory kinetic energy of Ti a b = 125 MeV. 
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FIG. 7: (Color online) pp and pn scattering observables versus the centre-of-mass angles 6 at an 
incident laboratory kinetic energy of Ti a b = 200 MeV. 
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